Semileptonic decays of double heavy baryons 
in a relativistic constituent three— quark model 



Amand Faessler 1 , Thomas Gutsche 1 , Mikhail A. Ivanov 2 , Jiirgen G. Korner 3 , Valery E. Lyubovitskij 

1 Institut fur Theoretische Physik, Universitat Tubingen, 

Kepler Center for Astro and Particle Physics, 
Auf der Morgenstelle 14, D-72076 Tubingen, Germany 

2 Bogoliubov Laboratory of Theoretical Physics, 
Joint Institute for Nuclear Research, 
14 1980 Dubna, Russia 

3 Institut fur Physik, Johannes Gutenberg- Universitat, 
D-55099 Mainz, Germany 

(Dated: July 3, 2009) 

We study the semileptonic decays of double heavy baryons using a manifestly Lorentz covariant 
constituent three-quark model. We present complete results on transition form factors between 
double-heavy baryons for finite values of the heavy quark/baryon masses and in the heavy quark 
symmetry limit which is valid at and close to zero recoil. Decay rates are calculated and compared 
to each other in the full theory, keeping masses finite, and also in the heavy quark limit. 

PACS numbers: 12.39.Ki, 13.30.Ce, 14.20.Lq, 14.20.Mr 

Keywords: relativistic quark model, double heavy baryons, semileptonic decays, form factors, decay widths 



1* 



* On leave of absence from Department of Physics, Tomsk State University, 634050 Tomsk, Russia 



2 



I. INTRODUCTION 

The semileptonic decays of double-heavy baryons provide yet another opportunity to measure the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element V c (,- This is particularly so since the transition matrix elements between 
double-heavy baryons obey spin symmetry relations in the heavy quark limit in addition to a model independent zero 
recoil normalization of the relevant transition matrix elements. In this paper we study current-induced transitions 
between double-heavy baryons in a fully relativistic constituent three-quark model. In the heavy quark limit we 
recover the spin symmetry relations among the form factors valid at zero recoil and close to zero recoil including their 
zero recoil normalization. Since the model is formulated in terms of finite values of the quark and baryon masses 
we are able to calculate the corrections to the spin symmetry relations and the zero recoil normalization valid in the 
heavy quark limit. 

Current-induced double-heavy baryon transitions have been analyzed in a number of model approaches. These 
include effective field theories based on heavy quark spin symmetry [l|, 0, H, HI , three-quark models [1, H 1 quark- 
diquark models [1, [9[ , and nonrelativistic QCD sum rules 0, El- The progress achieved up to now can be summarized 
as follows. In the heavy quark limit (HQL) or in the limit of infinitely heavy quarks (m c & — * oo) double-heavy baryons 
can be viewed as heavy-light meson-like states containing a heavy diquark in the 3 color state and a light quark in 
the 3 color state 0. In the HQL the spins of the light quark and heavy diquark system decouple. This gives rise to 
relations between different transition form factors involving double-heavy baryons in the quark-diquark picture [l[ 
as well as in the three-quark picture Q. In particular, working in the near zero recoil limit one can express all weak 
transition form factors between double heavy baryons through a single universal function r](oj) which depends on the 
kinematical parameter u> = v • v' where v and v 1 stand for the four-velocities of the initial and final double-heavy 
baryon, respectively. 

In Ref. 0] we have analyzed double heavy baryons for specific decay modes. We have restricted ourselves to spin 
1/2 to spin 1/2 transitions using the same relativistic constituent three-quark model (RTQM) [f| 12] as is being used 
in the present paper except that we now no longer have to rely on the impulse approximation. Differing from the 
approach of the present paper, in [f| we have treated the double-heavy baryons as bound states of a heavy 6-quark 
and a heavy-light (cq) diquark. In this paper we take double-heavy baryons to be bound states of a light quark and 
a double-heavy (be) diquark. In particular this means that the interpolating three-quark currents used in this paper 
have a different spin-flavor structure than the corresponding current in [5| . In the full theory this will lead to different 
predictions for the rates. We would like to emphasize, though, that, in the nonrelativistic limit, both currents are 
consistent with one another. 

The RTQM can be viewed as an effective quantum field theory approach based on an interaction Lagrangian of 
hadrons interacting with their constituent quarks. From such an approach one can derive universal and reliable 
predictions for exclusive processes involving both mesons composed of a quark and antiquark and baryons composed 
of three quarks. The coupling strength of a hadron H to its constituent quarks is determined by the compositeness 

r — —i I 1 ^ /2 

condition Zjj = [13l Il4j where Zjj is the wave function renormalization constant of the hadron. The quantity Z^ 
is the matrix element between a physical particle state and the corresponding bare state. The compositeness condition 
Zh = enables one to represent a bound state by introducing a hadronic field interacting with its constituents so that 
the renormalization factor is equal to zero. This does not mean that we can solve the QCD bound state equations but 
we are able to show that the condition Zh = provides an effective and self-consistent way to describe the coupling 
of a hadron to its constituents. One starts with an effective interaction Lagrangian written down in terms of quark 
and hadron variables. Then, by using Feynman rules, the ^-matrix elements describing hadron-hadron interactions 
are given in terms of a set of quark level Feynman diagrams. In particular, the compositeness condition enables one 
to avoid the problem of double counting of quark and hadronic degrees of freedom. The approach is self-consistent 
and all calculations of physical observables are straightforward. There is a small set of model parameters: the values 
of the constituent quark masses and the scale parameters that define the size of the distribution of the constituent 
quarks inside a given hadron. 

The main objective of the present paper is to present a comprehensive analysis of all possible current-induced spin 
transitions between double-heavy baryons containing both types of light quarks - nonstrange q = u,d and strange s. 
This involves the flavor transitions be — > cc and bb — > be where the transition be — > cc is treated as the generic process 
in the main text while the results for the transition bb — > be are mainly relegated to tables. The paper is structured as 
follows. First, in Sec. II we present interpolating three-quark currents with the appropriate quantum numbers of the 
double-heavy baryons. We then write down the corresponding Lagrangians defining the couplings of these currents to 
double-heavy baryons. Second, we briefly discuss the calculational techniques of how to calculate transition matrix 
elements generated by the Lagrangian functions. In Sec. Ill we consider the heavy quark limit of our transition matrix 
elements and recover the known heavy quark symmetry relations for the transition matrix elements between double 
heavy baryons as well as the appropriate zero recoil normalization of the form factors. In particular, we compare 
the results of the full finite mass calculation with the results derived in the HQL. Fourth, in Sec. IV we present our 
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numerical results which are compared to predictions of other theoretical approaches. In particular, we compare the 
results of the full finite mass calculation with the results derived in the HQL. Finally, in Sec. IV we present a short 
summary of our results. 

II. SEMILEPTONIC DECAYS OF DOUBLE-HEAVY BARYONS 

A. Lagrangian 

For the evaluation of the semileptonic decays we will consistently employ the relativistic constituent three-quark 
model. In the following we present details of the model which is based on an interaction Lagrangian describing the 
coupling between baryons and their constituent quarks. 

The coupling of a baryon -6(919293) to its constituent quarks 91, 92 and 93 is described by the Lagrangian 

£i%t( x ) = 9bB(x) JdxiJdx 2 Jdx 3 F B {x, xi, x 2 , x 3 ) J b (xi,x 2 ,x 3 ) + h.c. (1) 

where Jb(x\, x 2 , x 3 ) is the interpolating three-quark current with the quantum numbers of the relevant baryon B. 
One has 

Jb(xi,x 2 ,x 3 ) = e a ^T 1 q a 1 '{x 1 )q a 2 i {x 2 )CT 2 q^{x 3 ), (2) 

where the T 12 are sets of Dirac matrices, C is the charge conjugation matrix C = 7°7 2 and the (i=l,2,3) are color 
indices. Fb(x, x\,x 2 , x 3 ) is a nonlocal scalar vertex function which characterizes the finite size of the baryon. 

B. Vertex function 

The vertex function F B is related to the scalar part of the Bethe-Salpeter amplitude and characterizes the finite 
size of the baryon. To satisfy translational invariance the function Fb has to fulfill the identity 

F B (x + a,xi + a,x 2 + a,x 3 + a) = F B (x,xi,x 2 ,x 3 ) (3) 

for any given four-vector a . In the following we use a specific form for the vertex function 

F B (x,xi,x 2 ,x 3 ) = NsS^ix-^WiXi) <S> B [y^{x l - Xjf J (4) 

i=l i<j 

where $5 is a nonlocal correlation function involving the three constituent quarks with masses mi, m 2 , m 3 \ Nb = 9 
is a normalization factor. The variable uii is defined by Wi = mij \m\ + m 2 + m 3 ). The vertex function Q satisfies 
the translational identity ©■ 



The Fourier transform of the correlation function I J2( x i~ x j) ) can be calculated by using Jacobi coordinates. 

\i<3 / 

One has 

3 3 

$B(pi,P2,Ps) = N B J dxe-^f[ j dx^ iXi 8&(x - £ w iXi ) 4>s - x,)A 

i=l i=l i<j 

= (27r)U^{p~Y,P^B(-l 2 1 -l 2 2 ), (5) 

i=l 

where the Jacobi coordinates are defined by 

x\ = x + - ^ 6(2^2 + w 3 ) , 

X2 = x + 73 ^ lW3 + 7H &( 2wi + w ^ > 

x 3 = x - 4= £i(wi + w 2 ) + &(wi - w 2 ) . (6) 
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The corresponding Jacobi momenta read 

P = P1+P2+P3, 

h = 75^3(^1 +P2) ~ ^75 (m +w 2 )P3 , 

h = (2w 2 + w 3 )pi + (2wi + w 3 )p 2 + -jg (wi - w 2 )p3 , (7) 
where, according to Eq.((3]), WiXi — x. Since the function $>b[ \\^{ x i ~ Xj) 2 ) is invariant under translations, its 

i=l \i<j ) 

Fourier transform only depends on two four-momenta. The function l\ — 1%) in Eq. ([5]) will be modelled in our 

approach. The minus sign in the argument is chosen to emphasize that we are working in Minkowski space. A simple 
choice is the Gaussian form 

®b(-11 - ll) = exp(18(^ + Z 2 2 )/A|) (8) 

where the parameter Kb characterizes the size of the double heavy baryon. Since l\ and l\ turn into —If and —l\ in 
Euclidean space the form ([5]) has the appropriate fall-off behavior in the Euclidean region. 



C. Three— quark currents 

Double-heavy baryons are classified by the set of quantum numbers (J p ,Sd), where J p is the spin-parity of 
the baryon state and Sd is the spin of the heavy diquark. There are two types of heavy diquarks - those with 
Sd — (antisymmetric spin configuration [Q1Q2}) and those with Sd — 1 (symmetric spin configuration {Qi^})- 
Accordingly there are two J p = 1/2+ double-heavy baryon states. We follow the standard convention and attach 
a prime to the Sd — states whereas the £<j = 1 states are unprimed. Note that the J p = 3/2+ states are in the 
symmetric heavy quark spin configuration. In Table 1 we list the quantum numbers of the double-heavy baryons 
including their mass spectrum as calculated in @. 

We pause for a moment to discuss some of the features of the mass spectrum obtained in Q which are relevant for 
our calculation. One notes that there is a mass inversion in the (l/2 + ) mixed flavor states (Sf, c , Sj c ) and (0& c , f^ c ) in 
that M(E' bc ) > M(Eb c ) and M(n' bc ) > A/(fif, c ) contrary to naive expectation even though the heavy triplet diquark 
state has a higher mass than the heavy singlet diquark state in the model of i.e. m(bc; Sd — 1) > m(bc; Sd — 0). 
The inverted mass hierarchy is at the origin of the prime notation mentioned above. We mention that the inverted 
mass hierarchy is a feature of all models that have attempted to calculate the mass spectrum of double-heavy 
baryons @, [^, [ill, EE EH E3] ■ In particular, the inverted mass hierarchy implies that one can only expect substantial 
flavor-changing branching ratios for the two lowest lying states Sf, c and f2b c whereas the rates of the higher lying 
states S' bc and S£ c , and fl' bc and f2J c will be dominated by flavor-preserving one-photon transitions to the lowest- lying 
states Sb c and f2b c . It will be interesting to analyze the strength of one-photon transitions between the Sd — and 
Sd = 1 double-heavy baryon states which are forbidden in the HQL since, in the HQL, the photon couples to the 
light quark only. For finite heavy quark masses one-photon transitions between the Sd — and Sd — 1 double-heavy 
baryon states will occur at a somewhat reduced rate which, however, very likely will still exceed the flavor-changing 
weak decay rates of these states. 

We construct the interpolating currents of the double heavy baryon B q q 1 q 2 in the form of a light quark q ai coupled 
to a heavy diquark d < Q i Q 2 , viz. 

J qQl Q 2 = ?Q 1 Q,f 1 dg lQ *\ d^=s a ^ (q?CT%\ Q2 Q<?) . (9) 

We shall only consider currents without derivatives. With this restriction one can construct three interpolating 
currents for the (^ + ,0) states - the pseudoscalar J p , scalar J s and axial J A currents 

e a ^ q ai (QTC lb Q a 2 ^ , (10a) 

£ aia 2 a s 7 5 q a ^Q^CQ^ , (10b) 

e ms 7V 1 (QTCTol^QT) ■ (10c) 



JP 

• J qQiQ-2 



J, 



qQlQl 
A 

qQ\Ql 
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For the {\ + , 1) states one has a vector J v and a tensor J T current 



-v 

9Q1Q2 



e a lQ2 a 3 7 « 7 5 gai (q^C 7q Q^) , (11a) 

j'qQiQz = l^ 1 ™ 3 ^l 5 1 ai (Q?Ca^Q a A. (lib) 



Finally, for the (| + , 1) states one has the vector and tensor currents and Jj 

J 9 W , = £aia2a3 I* 1 (QTC^Q?) , (12a) 
JqQiQz, fi = -*e aia2a3 7V 1 (QTCa^Qf) . (12b) 

Note that any double-heavy baryon current in the form of a heavy quark coupling to a heavy-light diquark can be 
transformed to a linear combination of the above form of currents using a Fierz transformation. 

In the heavy quark limit the scalar current JQ 1 Q 2q vanishes, while the other currents become degenerate in the 
following way: 

JqQiQ2 = J? Ql Q 2 = ? aia2a3 4>? (C^Q 2 )> (13a) 
= JqQiQi = z aia2a3 (r Q \a 2 ar Q \), (13b) 

jJq.q, = j 7q>q. = ? aia2a3 i>r (r Q \wr Q \), w 

where ipq,Qi,Q2 are the upper components of the Dirac quark spinors and the <7j are Pauli spin matrices. Excluding 
the scalar current Jq 1 q 2(] , which vanishes in the HQL, we remain with two currents for each of the double-heavy 
baryon states which, as shown above, become degenerate in the HQL. It is therefore reasonable to take only one 
of the interpolating currents each. Our choice is to take the simplest current from each pair - the pseudoscalar 
current for the (i + , 0) states and the vector currents for the {\ + , 1) and (| + , 1) states. Note that the HQL coincides 
with the nonrclativistic limit. In the nonrelativistic limit our double heavy baryon (DHB) currents have a one-to- 
one correspondence to the naive quark model baryon spin- flavor functions which are displayed in Table 2. Further 
details on the naive quark model and how to evaluate the semileptonic current-induced transition amplitudes in this 
framework can be found in Appendix A. 

We now shall give explicit expressions for the three-quark currents which are needed for the calculation of the 
bcq — > ccq semileptonic transition amplitudes: 

J bcq = lf c y^ c) , d% c) = (ft-CltV*) , (14a) 

J bcq = dg^ri?, d a { l c) =e a ^^T^Cb a ^ (14b) 

where J = J^7° and T = ^T^j . The corresponding ccg-currents are obtained by the obvious replacement b — > c. 
Note that Cri^C = —T^ T . Using a rather suggestive notation we specify the coupling content of the DHB currents 
for the Zq 1 q 2 baryons using the form =>- <g> T^ d \ We thus consider the currents 



S cc =^ 7 Q 7 5 ®7a, Z* CC ^I®Y, (15a) 
S 6c => 7 Q 7 5 ®7a, Hk=>J®7", Z' bc =>I®^. (15b) 

We use the same set of currents for the double-heavy flQ 1 Q 2 -type baryons replacing the light u, d-quark by a s-quark. 



D. Normalization 



As described e.g. in [5[ we need the derivative of the mass operator for the double-heavy baryon 2q x q 2 in order to 
evaluate the coupling constants gs = gs QlQ2 ■ The mass operator is given by a two-loop Feynman diagram and reads 



n SQ lQ » 



h 



N„ 



(16tt 2 ) 2 
x tr 



d 4 fci f d i k 2 -=-2 2 2 

— V B ( — L 1 - l 2 ) 



r^Sifa + w lP )r ^s 2 (k 2 - W2P)\ r^s 3 (k 2 - ki 

j=(k 1 -k 2 ), h = ~^(ki + k 2 ), 



w 3 p 



(16) 
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where q\ = b or c, q 2 = c , q 3 — q and Nb c — 6, Af cc = 12. The expression Si(fc) = {m Qi — jt — ie)^ 1 denotes the free 
fermion propagator for the constituent quark with mass m qi . Integration momenta have been shifted in such a way 
so as to remove the external momentum from the vertex function. We have assigned outgoing momenta to the first 
and third outgoing quarks and an in-going momentum for the second quark. Since the p-dependence of the mass 
operator resides entirely in the propagators it is not difficult to calculate the derivative of the mass operator needed 
for the normalization condition 



(17) 



The latter condition is known as a Ward identity which is equivalent to the compositeness condition Zh = [H, EH • 
One obtains 



N, 



9192 



(16tt 2 



d A k X f d 4 k 2 -2 2 2 

J5T V M n h) 



{wi tr^S^h + w lP )^Si(h + w lP )t^S 2 {k 2 - w 2 p) 



x T {q) S 3 {k 2 - fci + w 3 p)f 



(9) 



w 2 tr 



T^S^h + w lP )T {d) S 2 (k 2 - w 2P )^S 2 (k 2 - w 2P ) 



x r^5 3 (fc2-fci + w 3 p))r (9) 



w 3 tr 



T^Sx(h + w lP )f {d) S 2 {k 2 - w 2 p) 



x T^S 3 (k 2 -h + w 3P )^S 3 {k 2 -ki + w 3 p)T («)} , 
where the double-heavy SQ 1 Q 2 -baryon is taken to be on its mass-shell. 



(18) 



E. Matrix elements of semileptonic decays 



The generic matrix element describing the semileptonic transitions Sj, c — > S cc reads 

(S cc (p')| J d 4 xe^c(x)0»b(x)\Z bc (p)) = (2n)U^(p - p< - q)u 3cc (pi)A»(p,p')u 3b M 

where 



(19) 



A"(p,j/) = 12: 



(16tt 



2V2 



x tr 



r<?s 4 (ii + hp')o^s 1 (k 1 +p-± P ')T ( £ ) s 2 (k 2 - \ v >) rWs 3 (fc 2 - fc^f 



=-(?) 

6c 



(20) 



Here O' 1 = 7^(1 — 7 ) and w\ = mi/(m,4 : + m 2 + m 3 ) with i = 4, 2, 3. For the masses one has mi = to;,, m 2 = 777,4 = rn c , 
771 3 = m q . The Jacobi momenta uii and uj\ are chosen as 



h 
h 



1 



1 



fcl - ^2 + W 3 p 

ki + k 2 + (w 2 - wi)p + q 



(21) 



fci — fc 2 + w' 3 p' 



1 

V6 



fcl + &2 + (t£?2 ~~ w a)p' 



Note that the expressions for the normalization and the vertex given in Eq. (|T5]) and Eq. (|2U)) . resp., are exact in the 
sense that they are obtained directly from the Lagrangian Eq. ([2]) for an arbitrary translationally invariant vertex 
function Fg such as the one defined in Eq. The two-loop integrals in Eqs. (|18l20p are invariant under translations 
of the loop-variables ki — > ki + bi (i=l,2) where bi are arbitrary momentum four-vectors. We have assigned the 
loop momenta such that the heavy quark limit can easily be taken. Calculational techniques for the two-loop quark 
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integrals are given in some detail in our previous publication One uses Schwinger's parametrization to raise the 
denominator factors into exponential factors. The tensor integrals are dealt with by using differential representations 
of the numerator factors. After doing as many loop integrations analytically as possible one ends up with four- 
fold parameter integrations for the derivative of the normalization factor and the transition form factors which are 
evaluated numerically. 



III. HEAVY QUARK SPIN SYMMETRY 



A. Structure of weak transitions in the HQL 



In the HQL the spins of the double heavy diquark and the light quark in a double heavy baryon decouple. At zero 
recoil and close to zero recoil this leads to spin symmetry relations among transition form factors between double- 
heavy baryons and a zero recoil normalization for the form factors. The spin symmetry is exact at zero recoil and close 
to zero recoil where the near-zero recoil region is specified later on. Near zero recoil one can choose the momenta as 

n^ccV^ + r M where r is a small residual momentum in the sense that r 2 ~ 0(1) when 



M — mb c v^ and ^ — 



Pi 



Since the final baryon is on mass-shell one has v ■ r = —r /2m e 

\{v-v') 2 = 2m 2 cc {l-w) 



0(l/m c ). This imposes a restriction on 
O(l). From the last equation one obtains 



Pi 

m. 

the kinematical variable w — v ■ v' since r 2 = 

w ~ 1 + 0(l/m 2 ). This situation is different compared with the case of baryons with a single heavy quark. A baryon 
with a single heavy quark possesses both a spin and a flavor symmetry. In what follows we will work near zero recoil 
in the sense that we neglect terms of 0(v ■ r). In order to keep things simple we assume that mj c = mi, + m c and 
m cc = 2m c . 

Using the above assumptions the heavy mass propagators simplify in the HQL. One has 



§b(h + p- 



1+ •} 



i 



2 — k\X> — it ' 
1± ft 1 



5c(fcl± ^' } 2 Thv-ie 



(22a) 



(22b) 



Due to the simple form of the heavy quark propagators in the HQL the defining equation for the coupling constant 
Eq. tfTo) simplifies considerably. One obtains 



1 = N Ql 



Hi 



■Q\Q 2 



(16tt 2 ) 2 



tr 



T (d) !+ tf f(d) !' 



d 4 fc! fd A k 2 -2. . T^S 3 (k 2 - fci)f (9) 

— *s(yo) : 



{—k\v — ie) 2 (k2V — ie) 



(23) 



A similar simplification occurs for the transition operator A fi (p,p') in Eq. (|2"0]) . One obtains 



A^{v,v') 



12 



: 5s c 



tr 



(16tt 2 ) 2 
d^ki [ d A k 2 



r (d) 



r 1 ' 

2 bc 



i(yo) ^b{vt 



T^S 3 {k 2 -k x )T 



(9) 
be 



{—k\v — ie)(—kiv' — ie)(k 2 v' — ie) 



(24) 



where 



y r = y(r) = -\{k x - k 2 f - \{k x + k 2 - rf . (25) 

It is not difficult to see that once we put v = v' (r = 0) in Eq. (f24|) the two-loop integral and the trace factor 
in (|2"4"|) reduce to the corresponding factors in Eq. (|2"3"|) . As a result all semileptonic transition matrix elements can 
be expressed in terms of a universal function rj(uj) normalized to 1 at zero recoil with u) = 1 in full consistency with 
the heavy quark spin symmetry results derived in Q . 
Neglecting terms of 0(y ■ r) one finds 
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where 



A^(s bc ( v ) 

A"(SU») 
A^(s bc ( w )^ 

A^M 



> H cc («') 

* H cc («') 

* S cc (t/) 



A^K^S^V) 



17(10) 
J(w) 



j(i) ' 
d 4 /ti 



$5(2/0) "M^) 



| 7 M 7 5 17(10) 



g vv ' <n{w) 



m q + (k 2 - ki)v 



(— k\v — ie) 2 {k2V — ie)(m 2 — (fe — fci) 2 — ie) 



(26a) 
(26b) 
(26c) 
(26d) 
(26e) 
(26f) 

(27) 
(28) 



By keeping the corrections of 0(v ■ r) one can obtain explicit model dependent expressions for the corrections to the 
spin symmetry relations Eq. (|26[) . These corrections will not be listed in the present paper. 

As an extra bonus of our dynamical treatment the bb — > be transition matrix clement can be related to the be — ► cc 
transition in the HQL. One simply has to replace m cc by mbb in the functional form of the universal function 17(10) in 
([27]). i.e. in the functions $5(2/0) and &B(y r ) appearing in (|28|) . 



B. Calculation of the universal function r){w) 

It turns out that one can derive a closed-form expression for the universal Isgur-Wise (IW) function r](w) if one 
uses a Gaussian ansatz for the three-quark correlation function $b as has been done in Eq. ([8]) . We use the Laplace 
transformation 



* B {z) = J dsd> L B ( S )e- sz (29) 


and the integral representation 

cxp ( _ TT^ X ) = fi V £i / dtl J dt2 cxp ( ~ Sl ((tl + + ® ~ 82 ^ + ' (30) 



-00 —00 



In terms of the variable w one obtains 



00 00 00 00 



J(w) = ~- J J J A2 ( mg + ^A^J J dh J dh ® b{Zw) $b(zi) ' (31) 



-00 -00 

where 

3 1 
A = -+a 3 , ® B (z) = d& B {z)/dz, 



lw-l\ 2 2 2 / 2 (1 + a 3 )(a? +ai«2 + af) + aia 2 a 3 \ 
= ( ii + m cc \j — - — j +t 2 + - la 3 m q + — I, (32) 
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and where m q is the mass of the light quark in the DHB. The integral (|3Tj) can be evaluated in closed form using the 
Gaussian ansatz for the correlation function Eq. (J5J. One obtains a rather simple form for the universal rj(w) function 



9 

r)(w) = cxp - 3 (to - 1) -p 2 



(33) 



The dependence on the light quark masses has disappeared due to cancelation effects between the numerator J(w) 
and the denominator J(l). 

For the slope p 2 of rj{w] defined by r)(w) = 1 — ,o 2 (ii; — 1) + ... one obtains 



di](w) 



dw 



Q cc 



(34) 



As mentioned before, the HQL results for the 66 — > foc-transitions can be obtained by the replacement m cc — > rribb 
in the IW function. Accordingly the slope of the IW function for the bb — > be transitions is obtained from (|34|) by 
the replacement m cc — ► rribb, i.e. the slope increases by the factor rn 2 b /m 2 c when going from the be —> cc case to 
the bb — > 6c transition if one uses the same size parameter A^ in both cases. One should stress that there exists a 
spin-flavour symmetry at zero recoil w = 1 giving 77(1) = 1, which means that the be — > cc transition is identical 
to the bb — ► 6c one. Close to zero recoil there exists only spin symmetry, because the IW functions for 6c — * cc and 
66 — > 6c transitions explicitly contain the flavor factors m cc and rribb, respectively. 



IV. RESULTS 



We now proceed to present our numerical results. We first present results on the semileptonic rates using finite 
heavy quark masses, i.e. we do not take the HQL for the transition form factors. In the finite mass case we present 
results for both the electron/muon mode and the r-mode. We then present results using the zero recoil form factors 
in Eq. (|26[) and the universal function r](ui) of Eq. (f3"3")) , i.e. we extent the validity of the near zero recoil approach 
to the whole kinematic range < q 2 < (Mi — M2) 2 . Also we present estimates of the width using the nonrelativistic 
quark model, in which, as described before, the wave functions have the same spin-flavor structure as our relativistic 
current considered in the nonrelativistic limit. Note, in the naive quark model we drop the <? 2 -dependence of the 
corresponding form factors and use only their values at q 2 = 0. We choose the Gaussian form Eq. © for the correlation 
function of the double heavy baryons. Our results depend on the following set of parameters: the constituent quark 
masses and the size parameter A^. The parameters have been taken from a fit to the properties of light, single and 
double heavy baryons in previous analysises [12 ]. One has 

m u ( d ) m s m c m b Ab , . 

0.42 0.57 1.7 5.2 2.5 - 3.5 GeV ' [ ' 

All our analytical calculations have been done using the computer program FORM [l|| . For the numerical evaluation 
we have used FORTRAN. 

In Table 3 we present results on the q 2 = values of the transition form factors Fy and F± in the naive quark model. 
In Table 4 we compare our finite mass results with the theoretical approaches of Ref. 0] and Q for the electron/muon 
mode. As indicated in Tables 4 and 5 we allow the size parameter Ab to vary in the range 2.5 < Ab < 3.5 GeV. The 
variation of our rates listed in Tables 4 and 5 reflect via error bars the variation of the size parameter A b ■ Note that 
a smaller value of Ab gives smaller rates and vice versa. In Table 5 we present detailed results on semileptonic rates 
of double heavy baryons in the exact finite mass approach for the electron/muon and tau-lepton modes. These are 
compared to the corresponding results in the heavy quark spin symmetry limit and in the naive quark model. One can 
see that the naive quark model gives larger rates due to the omission of form factor effects that would result from the 
nonlocal structure of semileptonic transitions. The results in the heavy quark symmetry (HQS) limit are calculated 
using the IW function (f3"3"|) with m cc — 2m c — 3.4 GeV for the 6c — > cc transitions and rribb = 2m& = 10.4 GeV for 
the 66 — > 6c transitions. 

In order to check on the quality of the HQS limit in Figs.l and 2 we show plots of the w-dependence of the form 
factor FY (w)/V2 for the exact finite mass result and the HQS result using the IW function rj(u)). As examples we 
take the two transitions Sf, c — > S cc and 5^6 — > S(, c . From Table 5 and Figs.l and 2 one can see that HQS limit is 
better justified for 66 — > 6c transitions, while there is a larger difference between the finite mass results and the HQS 
results for the 6c — * cc mode as is evident in Fig. 1 and Table 5. 



10 



V. SUMMARY 



In this paper we have analyzed the semileptonic decay of double heavy baryons using a manifestly Lorentz covariant 
constituent quark model approach. Our main results are summarized as follows: 

- We have derived results for the matrix elements of the semileptonic decays of double heavy baryons for finite 
values of the heavy quark/baryon masses and for the HQS limit which is valid at and close to zero recoil; 

- We have presented a detailed numerical analysis of the decay rates for the two (e//i)- and T-lepton modes in the 
exact finite mass approach and in the HQS limit; 

- We have compared our results with the predictions of other theoretical approaches. 

We hope that the results of this paper can be used to extract the value of the CKM matrix element V c b from future 
experiments on the semileptonic decays of double heavy baryons. 
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APPENDIX A: NONRELATIVISTIC QUARK MODEL: SPIN-FLAVOR WAVE FUNCTIONS AND 
SEMILEPTONIC DECAY CONSTANTS OF DOUBLE HEAVY BARYONS 



In this Appendix we present results on the q 2 = values of the double heavy transition form factors -^(0) and 
Fi (0) in the nonrelativistic quark model. As emphasized before the nonrelativistic quark model is based on the spin- 
flavor wave functions which arise in the nonrelativistic limit of the relativistically covariant double heavy three-quark 
currents with quantum numbers J p = ^ + and | + . The corresponding quark model spin-flavor wave functions are 
given in Table 2, where we use the following notation for the antisymmetric \A and symmetric (xs, Xs) spin wave 
functions with S z = +1/2 (spin projection on z-axis): 

xa = y§ {t (twt)}> xs = y^jt (n + it) - 2 mj, x* s = y^jm + m + m}- (ai) 

We derive the expressions for the semileptonic decay constants -F'i'(O) and F-^(0) of double heavy baryons using the 
master formulas: 

3 3 

F^(Q) =< B'\J2lIbc} {l) \B > and F?(0) =< B'\ Y,Whc\®\B > , 

i—l i=l 

where 03 is the z-component of the Pauli spin matrix and Ib c is the flavor matrix responsible for the b — ► c semileptonic 
transitions. In Table 3 we list the results for the SQ 1 Q 2 -type baryon decay modes. Corresponding results for the 
&,Q 1 Q 2 -type baryon decay modes can be obtained by replacing the light nonstrange quark q by the strange quark s. 



APPENDIX B: SPIN-KINEMATICS OF SEMILEPTONIC DECAYS 



In this Appendix we first write down covariant expressions for the current-induced transitions involving the (l/2 + ) 
and (3/2 + ) baryons which allows us to define sets of vector and axial vector invariant transition form factors. We then 
calculate all hclicity amplitudes expressed in terms of linear combinations of the invariant form factors. The advantage 
of using helicity amplitudes is that one obtains very compact expressions for the decay rates including lepton mass 
effects. We mention that the use of helicity amplitudes allows one also to derive very compact expressions for single 
angular decay distributions and for joint angular decay distributions of the decay products (see e.g. [2(| HH). 

The momenta and masses in the semileptonic decays of double-heavy baryons are denoted by 

Bi(pi, Mi) B 2 (p u M 2 ) + l(pi,mi) + Vl { Vvi 0) (Bl) 
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where p\ = p 2 + q and q = pi + p v . The matrix elements of the vector and axial vector currents J^ A ^ between the 
baryon states with spin 1/2 or 3/2 are written as 



Transition 



i + 



i+ . 

2 



m; = (b 2 \ j; \ Bl ) = u( P2 , S2 ) 7 „ W) i°» v ^rn (<?) + ^ (?) U ( P1 , Sl ) 



1" vV/„2\ 



Mi" 



M 1 



M A ^(B 2 \J A \B 1 )=u(p 2l s 2 )[ 1 ,F A (q 2 )- l a^F A (q 2 ) + ^ 



Transition ^ 



1+. 

2 ' 



Transition | 



i+ . 

2 ' 



Ml = (BZ\Jl\B 1 )=u"(p 2 , S2 )[g a ,Fr(q 2 )+ 1 ^F 2 V (q 2 ) 
+ P ^F?(q 2 ) + P ^F?(q 2 )} l5 u( Pl , Sl ) 

M 2 3 ^ ^ ~Mf 4 ^J u ^ 1 ' Sl ' ) 



= (S 2 |J^|Sr)=S( P2 , S 2)[^i / ('Z 2 )+7^^(9 2 ) 
+ ^W) + ^W)lw,»(p 1|fll ) 



M" 4 



= (B 2 |J^|BD = ufo, S2 ) [ft^FfV) + 7^^V) 

+ iT72 ^3 (9 ) + 7T2 ^4 19 ) u (Pi, s l) 



M 2 



M 2 



Transition § 



3 + 
2 



9a 9/3 / jpT// 2\ ■ 9f r y/ 2 \ . 9ju T^y/ 2x\ 

+ ga ^ Mi g ^ qa Fy(q 2 )y( PllS1 ) 

M A = (B^\J A \Bl)=u«(p 2 ,s 2 )[g ap (^F 1 A {q 2 )- l a^ 
f(7^( (Z 2 )-^f^(9 2 ) + ^(9 2 )) 

75U /3 0l,Sl) 



+ 



9a 9/ 3 

M? 



SV9£_-5^9a pi, 2 x' 



(B2a) 
(B2b) 



(B3a) 



(B3b) 



(B4a) 



(B4b) 



(B5a) 



(B5b) 



where u^ v = (i/2)(j fl j u — 7^7^) and all 7 matrices are defined as in Bjorken-Drcll. 

Next we express the vector and axial helicity amplitudes JJ^^ in terms of the invariant form factors F^' A ', 

where \w = t, ±1,0 and X 2 = ±1/2, ±3/2 are the helicity components of the W s shell and the daughter baryon, 

respectively. Since lepton mass effects are taken into account in this paper we need to retain the temporal component 
"t" of the four-currents J^' A - We need to calculate the expressions 



Hli w =Ml A (X 2 )e^(\ w ) 



(B6) 



Note that the helicity of the parent baryon (Ai) is fixed by the relation Ai = A 2 — Xw We shall work in the 
rest frame of the parent baryon B\ with the daughter baryon B 2 moving in the positive z-direction: p\ — (Mi , 0) , 
p 2 = (£ 2 ,0,0,|p 2 |) and q = (q Q , 0, 0, -|p 2 |), where E 2 = Q+/(2M 1 ), |p 2 | = y /^Ql/(2M 1 ), q = M 1 - E 2 and 
Q ± = (Mi ± M 2 ) 2 - q 2 . 
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The J = \ baryon spinors are given by 



«i(pi,±l) = V2Mil X± ) (B7) 



where x+ = ^ q J anc ^ = y 1 J two_com P onent P auu spinors. 
The J = I baryon spinors are defined as 



« M (p, a*) = ^<lA^|^*)e>, A)«(p, s) (B8) 

A.s 

where (lA^s||s*) is the projection matrix element of the spin | onto spin i; e^p, A) is the polarization vector and 
u(p,s) are the ususal J = \ spinors defined above. In particular, the J = | spinors with helicities A = ±3/2, ±1/2 
read: 

u "( p,± f) = c m(p» ±1 ) u (p» ± ^)i 

Mp>±<^) = "\/|^(^0) U (p,±i) +y|e AJ (p,±l) U (p,T^) (B9) 
The polarization vectors corresponding to the parent and daughter J = | baryons are given by: 

e(pi, 0) = (0, 0, 0, 1) , e(pi, ±1) = -^=(0, T l, -*, 0) , 

c*(P2,0) = -^(|p 2 |,0,0,£i) ) e*(p2,±l) = ^(0,Tl,i,0). (B10) 
The polarization vectors of the W ff-shcii are written as 

= 4f(^0,0,-| P2 |), 
e*"(±l) = -^=(0,±l,i,0), (Bll) 
e*"(0) = -^(|p2|,0,0,(to). 

vr 

They satisfy the conditions 



?„e*"(±l,0) = 0, g„e*"(t) = vV- (B12) 



Using the above formulas for the spin wave functions with definite helicities one can then calculate the helicity 
- H v — H A 

\ 2 A W A 2 \\v A 2 Aw' 



amplitudes H\ 2 \ w = HY , — i?^ 1 \ , where the vector and axial components for the different spin transitions are 



defined by: 



Transition 



1+ , i+. 



Hi 



HL = a\AF?M + + F± 



M 

M. 



r) 



2 1 2 X V 1 2 Ml 



where 



Hf n 



M± = M!±M 2 



4( 



1 + 3 Mi 



M_ - K 



2 2 \ 



M : 

M_ 



r) 



M1M2 



M 2 + M| - g 2 
2MiM 2 



/2MiM 2 (u; + l) 



^ Mr ) ' 



We 



g 2 = 2MiM 2 (uw 

M 2 + M| 
ax ~~ 2MiM 2 ' 



/2MiM 2 (w7- 1) 



aV, = 2^M 1 M 2 (t« - 1) , 

2 

= 2VM!M 2 (u; + 1) . 

2 

From parity or from explicit calculations one has for the \ 



1 + 



helicity amplitudes 



HI 



A9 , — Xv 



H 



v 



Hi 



■^2 j ~~ -^V 



-if 



A 2 ,A W ■ 
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1 + 3 + 

Transition t> —> § : 



itf t = -^aY t(w -l)(^M 1 -^M + + ^^(M lU ,-M 2 ) + J F < rA_) 
2 



ffjo = -y 3 aY (^(Miw - M 2 ) - F 2 y ( w + 1)M_ + F^(w 2 - 1)M 2 ) , 

Ht 2l = ±ai 2l (Ff-2F-( W -l)), 
From parity or from explicit calculations one has for the i + — > | + helicity amplitudes 



Transition | 



i+ . 



Hji = -Vfara--l)(^M 2 + ^M + -F 3 ^(M 1 -M 2W )-^), 
HY = -Jl*X (fY{M 1 - M 2 w) + F 2 v (w + 1)M_ - F 3 y (™ 2 - l)Mi) , 



F7 y — 1 ™^ P y 



fff t = /1 4 («> + 1) (^M 2 - itf M_ - ^^(Mi - M 2 w) - F^ ) 



2 

M 2 

2' 



(B16) 



(B18) 



*fi = ^4(-^ + 2F 2 >-l)), 
From parity or from explicit calculations one has for the | + — > i + helicity amplitudes 
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Transition | — > | : 

ttV _ ttV i-pV P V pV\ 

(B20) 

^ = ^Hf t (Fr,Fr,Fr) + l(w^-l)^Hf t (F 4 A ,F 5 A ,F e A ), 

= -f ^) + |(- 2 -l)^ 1 (^,^, M-^iw-l^Ff, 

Hi = -^HftF?,Ft,Ft)- ±^-1)^1?, 
HU = ^Hi{Ff,Ff,F£) - ±=ai{w-l)F 7 A , 

ttA _ tjA(tpA tpA t?A\ 
H \t - ~ H ^t\ t l 1*2 '^3 ) ! 

where 

2 

HY t (x,y,z) = aX t (xM- + z-^j , 

2 

HY (x,y lZ ) = aX (xM++yjj-}, 
HY 2l (x,y,z) = -aX 2l (x + y^±), 

(B21) 

2 

Hf (x,y,z) = ai (xM--yj-), 

ttA , N 4 / M_\ 

H hl {x,y,z) = a^-x + y—), 
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3 + 3 + 

From parity or from explicit calculations one has for the | — ► | hclicity amplitudes 



— A 2 ,-Aw A 2 ,A W ' 



A 2 , — A w 



= -H 



A 2 ,A W 



The decay width is given by the expression 



Mi 



r 

A Si — »s 



192tt 3 M? 



(B22) 



(B23) 



where N, 



Sl s 2 = 1 for si = 1/2 and 1/2 for s\ — 3/2; Tt Sl ^ S2 are the bilinear combinations of the helicity amplitudes: 
Hui = l^iil' + l^-i-il' + l^iol' + l^-iol 2 



n^s = |^ 1 | 2 + |^_i_ 1 | 2 + | J ff3 1 | 2 + | J ff_3_ 1 | 2 + |^ |0 | 2 + | J ff_i | 2 

+ ^|(3|^ t | 2 + 3|tf> t | 2 + | % | 2 + \H_^ + + \H_^ 

+ \H i0 \ 2 + \H_i f) 

Wi^i = I 2 + \H_1_1 I 2 + \H1_1 I 2 + I 2 + |^in| 2 + |-^-in| 2 

2 2 '2 2 '2 2 ' 2 J 2 J 
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+ ^(3|^ t | 2 + 3|ff_: t | 2 + \H kl \ 2 + + l^^l 2 + \H_ h ,\ 2 

+ l^iol 2 - 



H3 .3 
2 2 



l^ol 2 ) 

= l-ffii I 2 + \H_1_1 1 2 + l-ffsi I 2 + |.H"_3_ ] I 2 + \H1_1 1 2 + ii | 2 

2 2 '2 2 '2 2 

2 

+ \Hi \ 2 + \H_i \ 2 + \H i0 \ 2 + \H_, \ 2 + ^±(3\Hi t \ 2 + 3\H_ h 



3\H 3 J 2 + 3\H_ 



\H hl \ 2 + \H_ h _ 1 \ 2 + \Hs 1 \ 2 + \Hs ,\ 



+ \H, \ 2 + \H_s \ 2 + IH^I 2 + \H_ hl \ 2 + \H, \ 2 + \H_, a \ 2 ) 
In the zero recoil limit the expressions for the rates simplify considerably. One obtains 



r, 1 

2 2 



r, 3 

2 2 



= p w J dw ^w 2 - l(l+(w)F 2 (w) + l^(w)G 2 (w)^) , 
1 

^max 

= Pw J dw\/w 2 - 1 l 2 {w) G 2 {w), 



W max 

r^i = \p w [ dw^/w 2 - 1 l 3 (w) G 2 (w), 
22 2 J 
1 

U) mai 

l -Pw J dw^w 2 -l ^l+(w)F 2 (w) + l^{w)G 2 (w^j 



r 3 3 

22 



(B24a) 



(B24b) 



(B24c) 



(B24d) 

(B25a) 
(B25b) 
(B25c) 
(B25d) 
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where F = , G = Ff, r = M 2 /M u w = (l + r 2 )/(2r) and 



Pw ~ 12^ Ml r ' 
lf(w) = (wTl)(3w mM ±l-2 W ), 

o 

W - 



h{w) = 2(w + l)^w n 
l 3 (w) = 2(w + l)(w n 



6r 

(w 2 — l)r 

- w + jH" 

6 

i±H = ^ ^l = ( w )— y— ± (Wmax ± l)(w 2 - I ! j . ( B2(, . 
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Table 1. Classification and mass values of double-heavy baryons. 
Mass values are used from [j| except for the E cc mass which is taken from [l8l |. 



Notation 


Content 




s d 


Mass (GeV) 


1 — <cc 


n< ccX 


1/2+ 


1 


3.5189 






1/2+ 


1 


6.933 


"be 


af&cl 
y iy°j 


1/2+ 





6.963 


>— 'oh 


y i j 


1/2+ 


1 


10.202 


1 — 'CC 


y i^^j 


3/2+ 


1 


3.727 


"be 


y i j 


3/2+ 


1 


6.980 


'bfe 


y i j 


3/2+ 


1 


10.237 




s{cc} 


1/2+ 


1 


3.778 




s{6c} 


1/2+ 


1 


7.088 




s[6c] 


1/2+ 





7.116 


Qbb 


s{66} 


1/2+ 


1 


10.359 


^cc 


s{cc} 


3/2+ 


1 


3.872 




s{6c} 


3/2+ 


1 


7.130 


^66 


s{6fe} 


3/2+ 


1 


10.389 



Table 2. Double-heavy baryon wave functions. 



Baryon 


Wave function 


Baryon 


Wave function 


' — 'CC 


qcc xs 




sec xs 


^bb 


qbb xs 


Qbb 


sbb xs 




■±q(bc + cb) xs 


Qbc 


7js(fec + c6) xs 


"6c 


-^q{bc-cb) X A 


Q'bc 


-±^s{bc-cb) xa 


' — 'CC 


- qcc xs 


^cc 


- sec xs 


"66 


- qbb Xs 


Qbb 


- sbb x*s 


—6c 


- -^q(bc + cb) Xs 


QL 


- -^s{bc + cb) x*s 



Table 3. (0) and GY(0) in the nonrelativistic quark model. 



Quantity 


"i>c ' i — 'CC 


"6c * "cc 


"6c y "cc 


"6c 


* ' — 'CC 


"6c * "cc 


1 — '6c ^ ' — 'CC 




^bb — > "6c 


"66 — » "6c 


"66 ~~ * "6c 


"66 


~> "6c 


"66 — ► "6c 


"66 ~~ > "be 


*f(0) 


V2 














V2 


J\ A (o) 


i 

72 




2 
3 




2 


V2 
3 


s/2 
3 



Table 4. Semileptonic decay widths of double heavy baryons in units of 10 GeV. 
Comparison with other approaches in case of light leptons (e, /x) in the final state. 



Decay mode 


Ref. [9] 


Ref. [4] 


Our results 


Decay mode 


Ref. [9] 


Ref. [4] 


Our results 


— ~ 

—66 —6c 


1 r*^ 

l.OO 


1 Q9+0-25 


U.oU it U.oU 


S'66 — ' i'6c 


l.iU 


o i 4 +0.20 
Ai4 -0.02 


U.oO zn u.oz 


— 6c ' — cc 


Z.OU 


9 r 7 +0.26 

-0.03 


Z.1U I U. lU 


o, i o 

J "6c ' a "cc 
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Table 5. Detailed analysis of semileptonic decay widths of double heavy baryons in units of 10 GeV. 



Decay mode 


Exact results 


HQS limit 


NQM 


e, modes 


t mode 


Hfc6 — ► H6c 


0.80 ± 0.30 


0.46 ± 0.13 


1.33 ± 0.61 
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